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1. Three-loop fermion and gluon form factors 

In this Section we consider massless QCD and discuss the viitual coirections to the photon- 
quark and Higgs-boson-gluon vertex (see also Ref. |[T]] and references therein). It is convenient to 
decompose the vertex functions according to the Lorentz structure and define the form factors Fq 
and Fg via Tq = Y^Fq{q^) ,r^^ = (171 • q2 g^^ - q\ (72) Pgiq^) , where q = q\+q2 and qi (172) is the 
incoming (anti-)quark momentum in the case of Fq, and Fg depends on the gluon momenta q^ and 
q2 with polarization vectors £^{qi) and e^{q2)- Both for and Fg, which are obtained by applying 
appropriate projectors, we have = <?2 = 0- Some sample Feynman diagrams contributing to Fq 
and Fg are shown in Fig. |l|. 

For our calculation we have used two different setups. The basic idea of the first one has 
been described in Refs. |^ : integral representations for the coefficients of the master integrals aie 
derived. They depend on the exponents of the denominators of the integral under consideration 
and the space-time dimension d. In the recent years a procedure has been developed to evaluate 
the resulting parameter integrals in the limit of large d (see, e.g., Ref. [^]). Knowing sufficiently 
many expansion terms the coefficient function can be reconstructed since (for fixed exponents) it 
is a rational function in d. The evaluation of the three-loop vertex corrections profited quite a lot 
from the experience gained in the context of the evaluation of the four-loop two-point functions [Q] 
and the findings of Ref. In the latter paper it has been shown that the recurrence relations of 
«-loop three-point functions are equivalent to (n+ l)-loop two-point functions. 

The second method has only been applied to the singlet diagrams (see, e.g.. Fig. ^h)) con- 
tributing to Fq. It relies on the idea to combine the Laporta method with the Grobner bases 
technique which has been published in the computer code FIRE [||]. 

We parameterize the results for 7^^ and Fg in terms of the bare coupling which allows us to 
factorize all occurring logarithms of the form ln(2^/M^), where = —q^ > 0, and to cast the 
expressions in the form (x = q,g) 




(1.1) 



We refrain from listing the results in terms of general SU(3) colour factors, which can be found in 
Ref. [^], however, we present for illustration the finite' part of Fg in the case of QCD where it takes 




(a) (b) (c) 

Figure 1: Sample Feynman diagrams contributing to the Fg ((a) and (b)) and Fg (c) at three-loop order. 
Straight and curly lines denote quarks and gluons, respectively. 

'We refer to Refs. Ilffl for the divergent contribution. 
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the form 

^(3) 14423912 384479^2 370649^3 280069^4 49167^2^3 199263^5 
Fg = — ; 1 7— — i — 1 : — h 



243 108 18 32 4 lo 4 
4527765C6 „^ ^ / 11801309 42296^ 41018^ 
^56 ''^'''^ + ^'''^'^ + "4 U58 81— + ^1 



1055C4 16982^5 \ 2/^ 2239573 4^2 3376C3 349C4\ 

8 45 n 4374 3 81 18 ' ^'^ 



The two constants Xg i and X9 2 take the numerical values X9 1 ?a 1428.9963678666183591 and 
X9 2 ~ 528.0583±0.0326 where Z9 1 is available analytically ||TI[] andX9 2 is known numerically [0, 



lip with the indicated precision. 

The new NNNLO results for the form factors constitute building blocks for a number of ap- 
plications. Among them are the virtual corrections to Higgs boson production in gluon fusion, the 
Drell-Yan process and the two-jet cross section in e^e^ collisions. 

2. Three-loop matching of the vector current 

In contrast to the previous Section we consider here QCD with «/, = 1 massive and nj massless 
quarks and evaluate the matching coefficient between QCD and non-relativistic QCD (NRQCD) for 
the vector current. This quantity is important for phenomena where two heavy quarks are produced 
in electron-positron annihilation or a bound state of two heavy quarks decays into a lepton pair. 

The vector currents in QCD and NRQCD are given by jv = Qy^Q, and = (^'^ O^Xi where 
Q denotes a generic heavy quark with mass mg and ^ and % are two-component Pauli spinors for 
quark and anti-quark, respectively, and a*^ {k= 1,2,3) are the Pauli matrices. The two currents jv 
and f- can be used to compute vertex corrections with two on-shell quarks and momenta q\ and ^2 
(F,. and F,,). From the requirement that the results agree up to power-corrections in iuq defines the 
matching coefficient Cy 

Z2ry = c„Z2z;% + ..., (2.1) 



where Z2 denotes the on-shell wave function renormalization constant [12, O] and quantities with 



a tilde are defined within NRQCD. The ellipses in Eq. (2.1) represent terms of order 1/mg which 
are neglected. 

For the evaluation of Cy it is convenient to consider = {qi +^2)^ ~ 4mg and apply the 
so-called threshold expansion [14] to F,, which identifies the hard, soft, potential and ultra-soft 
integration regions. The latter three contributions are present both on the left- and right-hand side 
of Eq. (2J) and thus cancel out. Only the hard contribution where q^ = 4mg and which is only 



present in Fy has to be evaluated. This reasoning is based on the use of Dimensional Regularization 
which is crucial for the evaluation of higher order corrections to Cy since scaleless integrals are 
automatically set to zero. This concerns in particular Z2 and F,, which are both identical to one. 
The two-loop corrections to c,, have been evaluated in Refs. [15, [T^, 17] and the fermionic 



three-loop contribution in Refs. [18, 19 1. The setup used in Ref. \u% is completely automated 



It is based on a chain of programs which work hand-in-hand. The starting point is QGRAF [20] 
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which generates the ampUtudes for each Feynman diagram. In a next step q2e and exp [21] are 



used in order to identify the topologies and generate input expression in FORM [ ]22[ ] format. The 
same input file containing the description of the topologies needed by exp is also used in order to 
provide the input for crusher [ p3| ] which performs the reduction of all occurring integrals to a 
basic set, so-called the master integrals. The topology file is also used for providing the necessary 
input for FIESTA [|^, 25] which is employed for the numerical evaluation of all master integrals. 



There are several checks on the correctness of our results. As far as the analytical part is 
concerned the most important one is the evaluation of the Feynman diagrams for general gauge 
parameter t, : we have checked that the linear ^ -term drops out in our three-loop result. A strong 
check of the numerical part of our calculation is the change of the master integral basis, which is 
achieved analytically with the help of the integration-by-parts relations generated by crusher. 
In the new expression, which is a completely different linear combination of master integrals, the 
numerical evaluation is again performed with the help of FIESTA. 

In numerical form the result for c,, is given by (for /i = mg): 



Cv = 1-2.67— + (-44.55 + 0.41«/) ( — I 

+ (cv,g - 0.93(8)?!/, - 0.09?!/,?!/ + 120.75(1)?!/ - 0.82?i^) + &{a'^) . (2.2) 

Once the purely gluonic contribution Cv,^ is available, which is expected to be numerically dom- 
inant, Cy can be used in the analysis of the third-order cross section o{e^e^ — )• tt +X) close to 
threshold or the extraction of the bottom quark mass from T sum rules. 

3. Virtual NNLO corrections to Higgs production in gluon fusion 

In this Section we go beyond the effective theory which has been used in Section [T] in order 
to define Fg and consider the gluon-Higgs vertex with finite top quark mass. This constitutes an 
important contribution to the gluon fusion process which has the largest production cross section for 



Higgs bosons both at Tevatron and LHC. Whereas the NLO corrections are exactly known [26, 27, 



28] at NNLO until recently only approximations for infinitely heavy top quark have been available. 



Recently this gap has been closed and in the works [29, 3^, ^] the NNLO corrections to the cross 
section o{pp ^ H + X) incorporating the top quark mass dependence have been evaluated. The 



virtual corrections have been evaluated before in Refs. [^, In the following we describe in 



more detail the computation of Ref. [32] 



The virtual contribution to the partonic cross section can be cast in the form 

«' = «-(i + f a'" + (f)'«'^' + ...), (3.1) 

where the LO cross section is given by aLo = Gp ajfo{p,e) 5(1 — .x)/(288-v/2?r(l — £)), with 
X = Mfj /§. \/f is the partonic center-of-mass energy. The function /o and the analytical results of 
the first five terms in the p = Mjj/M} — > expansion for 5'^) and 5'^^ can be found in Ref. [^]. 
Sample diagrams contributing to 5^^^ are shown in Fig. ^(a). 

The setup is similar to the one used in Section ^ A difference is connected to the asymptotic 
expansion in the Umit <^ M} which is performed with the program e xp [ 2 1 ] and independently 



4 



Three-loop vertex corrections 



Matthias Steinhauser 




(a) 




(b) 



Figure 2: (a) Sample diagrams contributing to the NNLO virtual corrections to gg — > H. Solid lines 
represent the top quark with mass and curly (dotted) lines massless gluons (the Higgs boson), (b) Finite 
part of 5^^' as a function of p. The longer-dashed lines include successively higher orders in p. 



with an in-house Perl program []34|]. In this way the three-loop vertex integrals are reduced to 



one-, two-, and three-loop vacuum integrals which are treated with MAT AD [|35[], and to massless 
one- and two-loop vertex contributions which are reduced to master integrals (see, e.g., Ref. [^]) 



with the help of FIRE [37] and an independent program based of the Laprota-method. 

Explicit results for 5'-^^ and 5^^^ can be found in Ref. [32]. In this contribution we discuss the 
convergence properties in Fig. ^b) where the finite part of 5^^^ is shown as a function of p. One 
observes good convergence up to p 3 which corresponds to Me ~ 1 -VM, ?a 300 GeV. 
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